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Abstract: We study a number of different ingredients related to 9 dependence, the non- 
dispersive contribution in topological susceptibility with the "wrong" sign, topological sec- 
tors in gauge theories, and related subjects using a simple "deformed QCD". This model 
is a weakly coupled gauge theory, which however has all the relevant essential elements 
allowing us to study difficult and nontrivial questions which are known to be present in 
real strongly coupled QCD. Essentially we want to test the ideas related to the U{1)a 
problem in a theoretically controllable manner using the "deformed QCD" as a toy model. 
One can explicitly see microscopically how the crucial elements work. 
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1. Introduction and motivation 

In the present paper, we study the topological susceptibility in a "deformed" QCD, wherein 
we can work analytically. The goal is to understand a number of different ingredients 
related to 6 dependence, the non-dispersive contribution in topological susceptibility with 
the "wrong sign", topological sectors in gauge theories, and many related subjects using 
a simple "deformed QCD" . This model is a weakly coupled gauge theory, which however 
has all essential elements allowing us to study difficult and nontrivial questions which are 
known to be present in real strongly coupled QCD. Essentially we want to test the ideas 
related to the U {1)a problem formulated long ago [1, 2, 3], see also [4, 5, 6] in a theoretically 
controllable manner using the "deformed QCD" [7] as a toy model. One can explicitly see 
microscopically how all the crucial elements work. 

As a result of these studies we find that the QCD vacuum energy density contains a 
contribution, that is not related to the physical massive propagating degrees of freedom in 
the system, but rather originates from nontrivial topological sectors of the theory. These 
studies shed some light on microscopic mechanisms related to the resolution of the old 
U{1)a problem [1, 2, 3, 4, 5, 6]. While numerically the Monte Carlo simulations support 
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the picture developed in [1, 2, 3, 4, 5, 6], there were no analytical computations which 
would provide hints as to the microscopic nature of the contact term with the "wrong 
sign". Still, 30 years after the resolution of the U{1)a problem, the nature and the source 
of the contact term remains a mystery. The goal of this paper is to fill this gap, to identify 
the gauge configurations which saturate the contact term. 

The strategy for the present study is as follows. We start with analysis of the non- 
dispersive "contact" term in a completely solvable 2-dimensional gauge theory wherein 
there are not any physical propagating degrees of freedom. We explain how the contact 
terms emerges in this case. We then consider a simplified ("deformed") version of QCD 
which, on one hand, is a weakly coupled gauge theory wherein computations can be per- 
formed in theoretically controllable manner. On other hand, this deformation preserves 
all relevant elements such as confinement, degeneracy of topological sectors, nontrivial 9 
dependence, presence of non-dispersive contribution to topological susceptibility, and other 
crucial aspects, for this phenomenon to emerge. This "deformed" theory has recently been 
developed [7] and will be further discussed in section 3. Finally, we compute the contact 
term in four dimensional "deformed QCD" in section 4 where we interpret the contact term 
in terms of tunnelling events between different topological sectors of the gauge theory. We 
explain how the Ward Identities are automatically satisfied when light quark is inserted 
into the system. 

2. The contact term and degeneracy of topological sectors 

In this section we present an overview of the nature of the contact term which is not related 
to any physical degrees of freedom. We explain how we know about its mere existence and 
how the anomalous Ward Identities (WI) require its presence for the gauge theory to be 
self-consistent. We also give a simple two dimensional example explaining how this term 
emerges in gauge theories. Put differently, the nature of this "weird" contribution is entirely 
determined by the topological properties of the model rather than the physical propagating 
degrees of freedom of the system. Thus, we will find similar behaviour in theories with 
similar topological properties irrespective of the particularities of the theories. 

2.1 The contact term 

We start with definition of the topological susceptibility x which is the main ingredient of 
the resolution of the U{1)a problem in QCD [1, 2, 3], see also [4, 5, 6]. The necessity for 
the contact term in topological susceptibility x can be explained in few lines as follows. 
We define the topological susceptibility x in the standard way: 



where 9 is the conventional 9 parameter which enters the Lagrangian along with topological 
density operator q{x), see precise definitions below. The most important feature of the 
topological susceptibility X; for our present discussion, is that it does not vanish in spite of 
the fact that q = dfj_K^ is total divergence. Furthermore, any physical state of mass mc. 




(2.1) 
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momentum /c — )• and coupling (0|g|G) = cq contributes to the dispersive portion of the 
topological susceptibility with negative sign^ 



Xdispersive 





<0, (2.2) 



while the resolution of the U{1)a problem, which would provide a physical mass for the 
r]' meson, requires a positive sign for the topological susceptibility (2.3), see the original 
reference [3] for a thorough discussion. 



Therefore, there must be a contact contribution to Xi which is not related to any propa- 
gating physical degrees of freedom, and it must have the "wrong" sign, by which we mean 
opposite to any term originating from physical propagators, in order to saturate the topo- 
logical susceptibility (2.3). In the framework [1] the contact term with the "wrong" sign 
has been simply postulated, while in refs.[2, 3] the Veneziano ghost had been introduced 
to saturate the required property (2.3). This Veneziano ghost field is simply an unphysical 
degree of freedom with the "wrong" sign in the kinetic term such that it generates the same 
contact term when integrated out. It should be emphasized that these two descriptions are 
equivalent and simply two separate ways of describing the same physics and that in these 
two pictures, the claim that the "contact" term does not come from physical propagating 
degrees of freedom is manifest. 

It should be mentioned here that the "wrong" sign in topological susceptibility (2.3) is 
not the only manifestation of this "weird" unphysical degree of freedom. In fact, one can 
argue that the well known mismatch between Bekenstein-Hawking entropy and the entropy 
of entanglement for gauge fields is due to the same gauge configurations which saturate the 
contact term in the topological susceptibility in QCD [8]. In both cases the extra term with 
a "wrong" sign is due to distinct topological sectors in gauge theories. This extra term has 
non-dispersive nature, can not be restored from the conventional spectral function through 
dispersion relations, and can not be associated with any physical propagating degrees of 
freedom. 

2.2 Topological susceptibility and contact term in 2d QED 

The goal here is to give some insights on the nature of the contact term using a simple 
exactly solvable two dimensional QED2 [9]. We follow [8, 10] to discuss all essential 
elements related to the contact term. 

We start by considering two dimensional photodynamics (without fermions) which is 
naively a trivial theory as it does not have any physical propagating degrees of freedom. 
However, we shall argue that this (naively trivial) two dimensional photodynamics never- 
theless has a contact term which is related to the existence of different topological sectors 
in the theory. Thus, the presence of degenerate topological sectors in the system, which 

^We use the Euclidean notations where path integral computations are normally performed. 



Xnon— dispersive 



lim / d'^xe'''^T{q{x),q{0)}) >0. 



(2.3) 



-3- 



we call the "degeneracy" for short^, is the source for this contact term which is not related 
to any physical propagating degrees of freedom. 

The topological susceptibility x iii this model is defined as follows 



e2 



lim / d^xe'^"" {TE{x)E{<d)) , (2.4) 
where q = ^E is the topological charge density operator and 



^~ 47r2 



d^x q{x) = — d?x E{x) = k (2.5) 
27r J 

is the integer valued topological charge in the 2d U{1) gauge theory, E{x) = diA2 — (92^1 is 
the field strength. The expression for the topological susceptibility in 2d Schwinger QED 
model when the fermions are included into the system is known exactly [11] 



(2.6) 



e2 



d'{x)-^Ko{fi\x\) 



where fi"^ = e^/vr is the mass of the single physical state in this model, and KQ{fi\x\) is the 
modified Bessel function of order 0, which is the Green's function of this massive particle. 
The expression for x for pure photodynamics is given by (2.6) with coupling e = in the 
brackets^ which corresponds to the de-coupling from matter field ip, i.e. 



XE^M = ^,j d'x[5\x)]=^^. (2.7) 

The crucial observation here is as follows: any physical state contributes to x with negative 
sign 

Xdispersive ~ Imi _^2 _ ^2 < ^^'^l 

n " 

in accordance with the general formula (2.2) in four dimensions discussed previously. In 
particular, the term proportional —Kq{^\x\) with negative sign in eq. (2.6) results from 
the only physical field of mass /x. However, there is also a contact term J d^x [^^(x)] in 
(2.6) and (2.7) which contributes to the topological susceptibility x with the opposite sign, 
and which can not be identified according to (2.8) with any contribution from any physical 
asymptotic state. In the two-dimensional theory without a fermion (photodynamics), there 
are no asymptotic states since there are no possible polarization states, and so it is clear 



^Not to be confused with conventional term "degeneracy" when two or more physically distinct states are 
present in the system. In the context of this paper, the "degeneracy" references the existence of winding 
states \n) constructed as follows; T\n) = |n + 1). In this formula the operator T is the large gauge 
transformation operator which commutes with the Hamiltonian [T, H\ — 0, implying the "degeneracy" of 
the winding states |n). The physical vacuum state is unique and constructed as a superposition of |n) states. 
In path integral approach the presence of N different sectors in the system is reflected by summation over 
G Z in eqs. (2.12, 2.13, 2.14). 

the factor in front of (2.6) does not vanish in this limit as it is due to our definition (2.4) rather 
than result of dynamics 
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that the contact term (2.7) is not related to any physical propagating degree of freedom. 
Likewise, with a fermion included, there is one physical degree of freedom, yet we see also 
the additional "contact" contribution in (2.6). 

This term has fundamentally different, non-dispersive nature. In fact it is ultimately 
related to different topological sectors of the theory and the degeneracy of the ground 
state as we shortly review below. Without this contribution it would be impossible to 
satisfy the Ward Identity (WI) because the physical propagating degrees of freedom can 
only contribute with sign (— ) to the correlation function as (2.8) suggests, while the WI 
requires XQEoi'm = 0) = in the chiral limit m = 0. One can explicitly check that WI is 
indeed automatically satisfied only as a result of exact cancellation between conventional 
dispersive term with sign (— ) and non-dispersive term (2.7) with sign (-I-), 



XQED = ^ J 



5'(^)-|^i^o(^k|) 



[1 - 1] = 0. (2.9) 



e2 



4^2 



e2 1 



1 2 



e2 



47r2 



Therefore, contact term actually plays a crucial role in maintaining the consistency of the 
theory, because the WI can not be satisfied without it. While the exact formula (2.6) is 
known, it does not hint at the kind of physics responsible for the contact term with the 
"wrong sign", mainly what sort of field configurations should saturate the contact term. 
Below, we provide some insights on this matter. 

2.3 The origin of the contact term — summation over topological sectors 

The goal here is to demonstrate that the contact term in the exact formulae (2.6) and 
(2.7) is a result of the summation over different topological k sectors in the 2d pure U{1) 
gauge theory. The relevant "instanton-like" configurations are defined on a two dimensional 
Euclidean torus with total area V as follows [11], 

Ap = -§e,..'^, eEi^^ = '^, (2.10) 

such that the action of this classical configuration is 

1 /■ o o 27r2/c2 

2 J '^"^ = ^'-''^ 

This configuration corresponds to the topological charge k as defined by (2.5). The next 
step is to compute the topological susceptibility for the theory defined by the following 
partition function 

2 = ^ /pAe-^/'^'^^'. (2.12) 

All integrals in this partition function are gaussian and can be easily evaluated. The 
result is determined essentially by the classical configurations (2.10) and (2.11) since real 
propagating degrees of freedom are not present in the system of pure U{1) gauge field 
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theory in two dimensions. We are interested in computing x defined by eq. (2.4). In the 
path integral approach it can be represented as follows, 

= /^^ / d2^i?(^)^(0)e-^^'^'"^'. (2.13) 

This gaussian integral can be easily evaluated and the result is as follows [10, 8], 

^47r2fc2 2^ 

fcez 

In the large volume limit V ^ oo one can evaluate the sums entering (2.14) by replacing 
SfceZ ~^ f dk such that the leading term in eq. (2.14) takes the form 

XE^M = ^,-V.^^.^ = —^. (2.15) 

A few comments are in order. First, the obtained expression for the topological suscep- 
tibility (2.15) is finite in the limit V — )■ oo, coincides with the contact term from exact 
computations (2.6), (2.7) performed for the 2d Schwinger model, and has the "wrong" sign 
in comparison with any physical contributions (2.8). Second, the topological sectors with 
very large k ~ \/ e^V saturate the series (2.14). As one can see from the computations 
presented above, the final result (2.15) is sensitive to the boundaries, infrared regulariza- 
tion, and many other aspects which are normally ignored when a theory from the very 
beginning is formulated in infinite space with conventional assumption about trivial be- 
haviour at infinity. Last, but not least: the contribution (2.15) does not vanish in a trivial 
model with no propagating degrees of freedom present in the system. This term is entirely 
determined by the behaviour at the boundary, which is conveniently represented by the 
classical topological configurations (2.10) describing different topological sectors (2.5), and 
accounts for the degeneracy of the ground state. In this way, large distance physics enters 
despite the lack of physical long distance degrees of freedom. Furthermore, ee know that 
this term must be present in the theory when the dynamical quarks are introduced into 
the system. Indeed, it plays a crucial role in this case as it saturates the WI as (2.9) shows. 

We conclude this section by noting that the contact term in the framework [9] can be 
computed in terms of the Kogut-Susskind ghost by replacing the standard path integral 
procedure of summation over different topological sectors above as follows. The topological 
density q= m.2d QED is given by = (^)^ — '^4>i^ where (j) is the physical 
massive field of the model and 0i is the ghost [9]. The relevant correlation function in 
coordinate space which enters the expression for the topological susceptibility (2.4) can be 
explicitly computed using the ghost as follows 



1 1 

+ 



e \2 
2^7 



5\x) - ^Ko{fi\x\) 



(2.16) 
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where we used the known expressions for the Green's functions. The obtained expression 
precisely reproduces the exact result (2.6) as claimed. In the limit e — )• when the fermion 
matter field decouples from gauge degrees of freedom we reproduce the contact term (2.7, 
2.15) which was previously derived as a result of summation over different topological 
sectors of the theory. The non-dispersive contribution manifests itself in this description in 
terms of the unphysical ghost scalar field which provides the required "wrong" sign for the 
contact term. These two different descriptions are analogous to the same two computations 
in four dimensions mentioned in section 2.1, with and without the Veneziano ghost, and 
again we emphasize the equivalence of the two. In the picture wherein the contact term 
is saturated by a ghost field we see again how the contact term is related not to physical 
propagating degrees of freedom. 



3. Deformed QCD 

Next we discuss the "center-stablized" deformed Yang-Mills developed in [7] and references 
therein, before moving on to a discussion of the topological properties of this theory in 
section 4. In the deformed theory an extra term is put into the Lagrangian in order to 
prevent the center symmetry breaking that characterizes the QCD phase transition between 
"confined" hadronic matter and "deconfined" quark-gluon plasma. Thus we have a theory 
which remains confined at high temperature in a weak coupling regime, and for which it is 
claimed [7] that there does not exist an order parameter to differentiate the low temperature 
(non-abelian) confined regime from the high temperature (abelian) confined regime. We 
follow [7] in deriving the relevant parts of the theory. 

3.1 Formulation of the theory 

We start with pure Yang-Mills (gluodynamics) with gauge group SU (N) on the manifold 
X with the standard action 

S^^ = [ d'x -l^tr {FlAx)\ , (3.1) 
and add to it a deformation action, 

AS= f d^x -^P[Jl(x)], (3.2) 
built out of the Wilson loop (Polyakov loop) wrapping the compact dimension, 

^^(x) =P p^"^^'* ^4(^-^4)] . (3.3) 

The "double-trace" deformation potential -P [fi] respects the symmetries of the original 
theory and is built to stabilize the phase with unbroken center symmetry. It is defined by 

[N/2i 

P[n]^ a„|tr[0-]p. (3.4) 

n=l 
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Here YN/2\ denotes the integer part of N/2 and {a„} is a set of suitably large positive 
coefficients. 

The first term of P [17], proportional to |tr (with a sufficiently large positive coef- 
ficient), will prevent breaking of the center symmetry from Z^v to Zi with order parameter 
(tr but will not prevent tr [fi-^j from developing a vacuum expectation value so that 
it will not prevent the center symmetry breaking from Z^v to Z2 with order parameter 
(tr[0^]). The term proportional to |tr [O^] | however does prevent such a symmetry 
braking. Likewise, for each other subset Zp of Z^r (with mod p = 0), there needs to 
be a corresponding term, proportional to |tr[$7P]|^, in the deformation potential. This is 
the reason for including terms up to |tr [fil-^/^-l] | . Note that for real life QCD the gauge 
group is SU (3) and so only one term is necessary, 

PP=a|trP|^ (3.5) 

In undeformed pure gluodynamics the effective potential for the Wilson loop is mini- 
mized for Vl an element of Zjy. The deformation potential (3.4) with sufficiently large {a„} 
however changes the effective potential for the Wilson line so that it is minimized instead 
by configurations in which tr [0"] = 0, which in turn implies that the eigenvalues of O 
are uniformly distributed around the unit circle. Thus, the set of eigenvalues is invariant 
under the Zjv transformations, which multiply each eigenvalue by e^'^'^^l^ (rotate the unit 
circle by k/N). The center symmetry is then unbroken by construction. The coefficients, 
{a„}, can be suitably chosen such that the deformation potential, -P[ri], forces unbroken 
symmetry at any compactification scales [7], but for our purposes we are only interested 
in small compactifications (L <^ where L is the length of the compactified dimension 
and A is the QCD scale). At small compactification, the gauge coupling at the compacti- 
fication scale is small so that we can work in a perturbative regime and explicitly evaluate 
the potential for the Wilson loop due to Quantum fluctuations as in [12, 7]. The one-loop 
potential is 

V[^= I d^x^V^x)], (3.6) 

with 

7 °° 1 

n=l 

The undeformed potential for the Wilson loop (3.7) is minimized when $7 is an element 
of the center, Z^r, so that Q. = g^Trjfc/A^^ The deformation potential (3.4) must therefore 
overcome this one- loop potential and force O to not choose one particular element of Z^v. 
We must choose the coefficients {a^} to be larger than 2/ (vr^n^). A simple choice is 
an = 4/ (vr^n^). With this choice, the full one-loop effective potential for the Wilson loop 
is minimized for tr [f]"'] = for all n 7^ mod N , indicating unbroken center symmetry. 

3.2 Infrared description 

As mentioned in the previous section, we are interested in the regime in which the compact- 
ification size is small, L <^ A~^, and so the gauge coupling is small a the compactification 
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scale, (1/-^) 1- So, in our deformed theory, the combined effective potential for the 
Wilson loop is the sum of (3.7) and (3.4), which is minimized by field configurations with 



VI = Diag (l, e^-V^, e^^^^/^, . . . , e2-{^-i)/^) , (3.8) 

up to conjugation by an arbitrary element of SU{N). The configuration (3.8) can be 
thought of as braking the gauge symmetry down to its maximal Abelian subgroup, SU (N) — >• 
U (1)^~^. In the gauge in which is diagonal, the modes of the diagonal components of the 
gauge field with zero momentum along the compactified dimension describe the U{1)^~^ 
photons. Modes of the diagonal gauge field with non-zero momentum in the compactified 
dimension form a Kaluza-Klein tower and receive masses that are integer multiples of 27r/L. 
The remaining off-diagonal components of the gauge field form a Kaluza-Klien tower of 
charged VF-bosons which receive masses that are integer multiples of 2tt/NL. Then the 
lightest VF-boson mass, myy = 2tt/NL, describes the scale below which the dynamics are 
effectively Abelian. 

As described in [7] , the proper infrared description of the theory is a dilute gas of N 
types of monopoles, characterized by their magnetic charges, which are proportional to 
the simple roots and affine root of the Lie algebra for the gauge group U{1)^. Although 
the symmetry breaking is SU{N) — > U{1)^~^, it is simpler to work with U{1)^ and, as 
we will see, the extra degree of freedom will completely decouple from the dynamics. The 
extended root system is given by the simple roots, 

ai = (1, -1,0, . . . ,0) = ei - 62, 

02 = (0, 1, -1, ... ,0) = 62 - 63, 



aN-l = (0, . . . , 0, 1, -1) = CN-l - CAT, 



and the affine root. 



aN = (-1,0, . . . ,0, 1) = eTv - 61. 
We denote this root system by Aaff and note that the roots obey the inner product relation 

For a fundamental monopole with magnetic charge aa € Aafr, the topological charge 
is given by 



Q 



,x5i 167r^ 



and the Yang-Mills action is given by 



Sym 



1 



^ 1 



>xS 



(fx -^tr 
1 2^^ 



1 



(3.10) 



\Q\. (3.11) 



The second equivalence hold because the classical monopole solutions are self dual, [12] 



Ffxv — F^ij. 
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For an antimonopole with magnetic charge —aa, the Yang-Mills action is the same (3.11) 
and the topological charge changes sign, Q = —1/N. 

So the infrared description, at distances larger than the compactification length L, 
is given by a three dimensional dilute monopole gas with N types of monopoles (and so 
N types of anti-monopoles) interacting by a species dependent Coulomb potential with 
interactions defined by the inner product (3.9), 



Va,bir) = L 



27tY {±aa) ■ {±ai,) 
g J 47r|r| 



27ry 2^ 



b-i 



Ja,b+1 



9 J 



Att Irl 



(3.12) 



where the overall sign is plus for a monopole-monopole or antimonopole-antimonopole 
interaction and minus for a monopole-antimonopole interaction. For a given monopole 
configuration with n^'^^ monopoles and n^""^ antimonopoles of types a = 1, . . . , A^, at posi- 
tions r[.°\ k = 1, ... ,71°" and f^'^^ k = 1, . . . , n" respectively, the three dimensional C/(l)^ 
magnetic field is given by 



N 



„ , , \ ^ 27r 

a=l ^ 



k=i 47r X — r 



(a) 



-(a) 



1=1 47r X — r 



(a) 
I 



Letting 



(a) 



for k < n^'^^ 
for k > n^'*) 



+1 for k < 
-1 for A; > 



(3.13) 



(3.14) 



we can write (3.13) in a more compact form, 



a=l ^ 



eq 

k=l 



(a) 
k 



47r 



X 



(3.15) 



The action for such a monopole configuration is a combination of the monopole self-energies 
and the Coulomb interaction potential energies for each pair of monopoles. 



where 



N 



5MG = SselfE^^^"^+'5mt> 



a=l 



(3.16) 



2tt^L 



N 



Jint 



E "a ■ «b 



a,b=l 



A: 



A:=l /=1 



(3.17) 



and 



G(r) 



47r Irl 



(3.18) 
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The canonical partition function is then given, as usual, by a sum over all possible monopole 
configurations with a statistical weight , 



Z 



N 



\{d^l 

a=l 



{a) g-^in 



(3.19) 



with measure 



k=l 1=1 

The monopole fugacity, C, describes the density of monopoles and is given by, 



(3.21) 



where the C factor is the one-loop functional determinant in the monopole background as 
described in the appendix of [7]. 

Next we show explicitly that the above monopole partition function (3.19) is equivalent 
to a Sine-gordon partition function which describes the proper ^-dependence for the QCD 
vacuum. 

3.3 Monopole sine-Gordon equivalence 

The sine-Gordon partition function for this model describes a three dimensional A^-component 
real scalar field theory, given by 

N 



with 



dual 



z 



-■S'dual ["■] 



a=l 



Considering the cosine term, 



exp 



■ N 

c / d'x y . 



2L \2tt 



N 



N 



{Vaf - C^cos(Qa • cr) 



a=l 



(3.22) 



(3.23) 



cos • a 



TT exp C / d^x cos(aa • a) 



N 

n 

a=l 



exp 



d^x (e*""-'^ + e"^""-'^) 



(3.24) 



we can apply the power series representation for the exponential, = '^x"'/n\, and get 



N 



(/./....) =n E 



(C/2) 



Mia) M(") 



n 



"=1 I M(<')=0 m=0 

We then make use of the binomial theorem 

n 

ix+yr = j2 



d Xyyi 



giaa-cr(xm) _|_ ^-iaa-(T{Xm) 



(3.25) 



k=0 



n 

k 



x" ''y'' with 



n 



k / {n-ky.kV 



(3.26) 
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and arrive at, 



N oo 



(/A..) = n E 

a=l Af{<i)=o 



E 

m=0 



(C/2) 



" TV 












exp 




.a=l 






a=0 k=0 











Si 


d^Xk 


1 exp 










fc=l 



(3.27) 

where d/i'-'^^ is given in (3.20). Thus, inserting (3.27) into the sine-Gordon partition function 
(3.22), we have 



Z = [Vaa dfi^'''> exp {-^ J , 



a=l k=l 



where 



/3 



1 / 9 



L \2tt. 

Treating the last term in the exponent as a source term. 



(3.28) 
(3.29) 



/3 



(3.30) 



a=l k=l 

and completing the square with the shift ct(x) — t- o"(x) + f d^y G(x — y) J(y), we have, 

N 



.a=0 



exp 



d^x / dS [J(x)G(x - y) J(y)] 



(3.31) 



in which is the functional determinant 



N 



.a=0 



exp 



d^x (Vo-)^ 



(3.32) 



The above determinant does not contain any of the relevant physics and is just a con- 
stant prefactor that will drop out of any calculation of operator expectation values in the 
monopole ensemble. Finally, inserting the expression for the source (3.30), the partition 
function becomes. 



Z = Z^ 





" N 




/ 




exp 




.a=0 





9 ah=\ k=\ «=0 



(3.33) 



which is the partition function for the monopole gas from (3.19). 
Next, including a ^-parameter in the Yang-Mills action. 



'S'yM — ^ ^YM + 



1 



ix5 



1 167r2 



tr 



(3.34) 
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with F^^ = e^'^'^'^ Fpcj, multiplies each monopole fugacity by e*^/^ and antimonopole fugac- 
ity by e"*^/^. In the dual sine-Gordon theory this inclusion is equivalent to shifting the 
cosine term so that^ 



The 6 parameter enters the effective Lagrangian (3.35) as 9/N which is the direct con- 
sequence of the fractional topological charges of the monopoles (3.10). Nevertheless, the 
theory is still 27r periodic. This 27r periodicity of the theory is restored not due to the 2n 
periodicity of Lagrangian (3.35). Rather, it is restored as a result of summation over all 
branches of the theory when the levels cross at = Tiirnod 27r) and one branch replaces 
another and becomes the lowest energy state. Indeed, the ground state energy density is 
determined by minimization of the effective potential (3.35) when summation X^^g^ over 
all branches is assumed in the definition of the canonical partition function (3.19). It is 
given by 



where V is 3d volume of the system. Eq. (3.36) shows that in the limit y — t- oo cusp 
singularities occur at the values aX. 9 = Tr{mod 2tt) where the lowest energy vacuum state 
switches from one analytic branch to another one. The first derivative of the vacuum 
energy, which is proportional to the topological density condensate, is two-valued at these 
points. This means that whenever 9 = ^{mod 27r) we stay with two degenerate vacua in 
the thermodynamic limit. If, on the other hand, the thermodynamic limit is performed 
for a fixed value of 9, any information on other states is completely lost in eq. (3.36). 
Correspondingly, the 27r periodicity in 9 is also lost in infinite volume formulae. We have 
no chance to know about additional states when we work in the infinite volume limit 
from the very beginning. As a result, usual V = oo formulae become blind to the very 
existence of a whole set of different vacua, which is just responsible for restoration of the 
2-K periodicity in 9. The model under consideration explicitly supports this pattern in 
deformed QCD where all computations are under complete theoretical control. 

Such a pattern is known to emerge in many four dimensional supersymmetric models, 
and also gluodynamics in the limit N = oo. It has been further argued [13] that the same 
pattern also emerges in four dimensional gluodynamics at any finite N . We follow, in fact, 
the technique from [13] to arrive at (3.36) in analyzing the 9 periodicity of the theory. The 
same pattern emerges in holographic description of QCD [14] at = oo as well. 

*We note in passing that there is a typo in [7] in sine-Gordon representation which is corrected here. 
Also, it has been stated (wrongly) in ref. [7] that the sine-Gordon Lagrangian is 2tt periodic as a result of a 
symmetry. This statement is incorrect, as the claimed symmetry is not in fact a symmetry of the theory, 
such that 9 parameter enters the Lagrangian as 6/N as it should. To check this insert cr = and notice 
that the 6-dependence is explicitly different after the transformation suggested in [7]. 




(3.36) 
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Finally, considering the expectation value 



±ia6-(T(y)\ 



1 

~z 

Z 



N 



.a=l 



exp <^ -/3 / d^x 



1 2 ^ 

- (Vcr) + J- <T=F -(5(x - y)afc • cr 



" N 




n ^/^^"^ 


g-SMG gxp 


.a=0 





47r^L (a) 



^ a=l fc=l 



(3.37) 

we note that the operator e^°'a-T{:x.) jg ^.j^g creation operator for a monopole of type a at x, 
i.e. 



(3.38) 



Likewise, the operator for an antimonopole is A^(j(x) = e *"<i cr(x)_ 'pj^g expectation values 
of these operators in fact determine the ground state of the theory. 

3.4 Mass gap 

The cosine potential in the sine-Gordon action (3.23) gives rise to a mass term for the dual 
scalar fields. Expanding the potential around the minimum <t = up to quadratic order 
and rescaling <t — )■ L{2tt)'^ / a to put the kinetic term into canonical form, gives (up to a 
constant term) 



1 ^ 

V{<T)^-mlY,{^a+l-CJi)\ 



(3.39) 



a=l 

with 

The above mass term is diagonalized by the discrete Fourier transform 

N 

'N 



(3.40) 



1 ^ 

/AT ^-^ 



(3.41) 



becoming 



a=0 



N 



V{ct) 



1 2 I ~ |2 

a=l 



(3.42) 



where ma = m„ sm^ira/N). So the only scalar field which remains massless is the A^th field, 
which is the field associated with the affine root. Inserting the discrete Fourier transform 
(3.41) into the full cosine potential however shows that the A^th field drops out of the 
cosine potential completely, so although it remains massless, it completely decouples from 
the theory and does not interact with the other components at all. 
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4. Topological susceptibility in the deformed QCD 



Next we consider the topological susceptibility in the deformed theory discussed in the pre- 
vious section in both the monopole and sine-Gordon formalisms. We define the topological 
density g(x) and topological charge Q by 



1 



d X iv 



and as in (2.1) the topological susceptibility x is given by 



X = L m-n / d^x e^^ "" (g(x)g(O)) 



L / d'x fe(x)] 



(4.1) 



(4.2) 



First, in next subsection we compute the topological susceptibility directly, using the 
monopole gas representation. As the next step, we reproduce our results using sine Gor- 
don representation of the theory. Finally, we compute the topological susceptibility when a 
single massless quark is introduced into the system. Essentially, the goal here is to discuss 
the same physics related to non-dispersive contact term, topological sectors and all that in 
deformed QCD, in close analogy to our discussions in 2d QED in section 2. 

4.1 Topological susceptibility in the monopole picture 

In order to compute the functional form of the topological susceptibility in the monopole 
theory we consider the topological density, 



IF^^F^'' 



-1 



Ajk4 



9 
47r2 



a=l 



N 

5:(4'^))[v.bW(x)], 

a=l 

where the U{\)^ magnetic field, = e^^^^Fjk/2g is given by 



V F^"' F 



(a) j^{a) 
j4 



(4.3) 



2tt 



B('^)(x) = —a, 



X — r 



E 

A;=i 47r |x — r 



(a) 



X — r 



(a) 



E 

k=i 47r |x — r 



(a) 



(a) 



(4.4) 



and (^aI^^^^ is just the expectation value of the diagonal gauge fields in the compact direc- 
tion, 



A 



(a) 



2-K 



(4.5) 



The above fi'^ are the fundamental weights for the SU (N) algebra and are defined by 

^-.a,^\5^al = 6l (4.6) 
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Inserting the magnetic field for the monopole ensemble into the topological density expres- 
sion (4.3) and applying Gauss's theorem to the result, we arrive at 



TV 



N 1 
^ Tn 

a=l [fe=l 1=1 

a=l k=l 



(4.7) 



which obviously gives the proper topological charge for a single monopole or antimonopole, 
Q = ±1/N. The topological density operator q(x) has dimension four as it should. 

The expectation value (g(x)g(O)) is the topological density operator (4.7) evaluated at 
each point inserted in the partition function (3.19), 

N 



(g(x)g(O)) = ^ [ Hci/x^'^) [g(x)<?(0)]e-^- 

a=l 

.AT N Af mC-) 

= H""'" E E E K°'<3!"*(r[°' -x)«(r!")] 

a=l a,b=l k=l 1=1 

~ZIW / ^/^EE[^™'3n 5(r„-x)(^(r„)]e-^-' 

m n 

j dn I 5(x) 5(r„)e"^-' + E E ['3™^" ^("^"i - x)(5(r„)] e" 



c 



{5(x)-0(C)}, 



iVL2 

(4.8) 

where we have condensed the indices to just m and n which run over each monopole in the 
ensemble. In the above expression, the double sum of delta functions gives a set of terms 
in which each pair of monopoles in the ensemble are moved to the points x and and 
computes the partition function given that arrangement. The monopole gas experiences 
Debye screening so that the field due to any static charge falls off exponentially with 
characteristic length . The number density J\f of monopoles is given by the monopole 
fugacity, ~ so that the average number of monopoles in a "Debye volume" is given by 



The last inequality holds since the monopole fugacity is exponentially suppressed, C ~ 
e~^/^ , and in fact we can view (4.9) as a constraint on the validity of our approximation. 
The statement here is that inserting or removing a particular monopole will not drastically 
affect the monopole ensemble as a result of condition (4.9), so that we can compute expec- 
tations of operators in the original ensemble without considering the back-reaction on the 
ensemble itself. The factor overall factor of C, and additional factor in C'(C) in formula 
(4.8) appears because each monopole we remove from the ensemble leaves a factor of i^/iV^") 
in the monopole measure, and there are A'^('^) such terms for each type of monopole so that 
we are left with just a factor of 



-16- 



The computed non-dispersive contribution (4.8) to the topological susceptibility in de- 
formed QCD has exactly the same structure we observed in two dimensional QED discussed 
in section 2. In particular, it is expressed in terms of a (5(x) function, and it has the "wrong 
sign" similar to (2.7). Furthermore, this contribution is not related to any physical prop- 
agating degrees of freedom, but rather, it is determined by degenerate topological sectors 
of the theory. The corresponding "degeneracy" is formulated in terms of monopoles which 
essentially describe the tunnelling transitions between those "degenerate" sectors, see sec- 
tion 4.3 for more comments on the physical meaning of the formula. If we neglect a small 
term 0{C) in formula (4.8) we arrive to the following final expression for the topological 
susceptibility in deformed QCD without quarks 

XYM = ^j d^x [S{^)] = ^. (4.10) 

It has dimension four as it should. This expression is a direct analog of eq. (2.7) derived 
for two dimensional QED. The same formula (4.10) could be computed much easily in the 
dual sine-Gordon theory by differentiation of the ground state energy density (3.36) with 
respect to the 9 as general expression (2.1) states 



Xym{0 = 0) 



89^ 



^, E^u.{9 = 0) = -^. (4.11) 



Agreement between the two computations can be considered as a consistency check of 
our approach in weakly coupled regime. One can explicitly see that the general relation 
Xym ~ -£'min(^ = 0)/A^^ holds for the deformed QCD as a result of 9/N dependence in eq. 
(3.35). The difference with real strongly coupled QCD is that the vacuum energy scales as 
iV^ in QCD rather than ~ iV in (4.11). 



4.2 Topological susceptibility in the presence of the light quarks 

Our goal here is to introduce a single massless quark ip into the system to see how the 
topological susceptibility changes in this case. We anticipate that the emerging structure 
should be very similar to (2.6) as the topological susceptibility must vanish in the presence 
of massless quark in the system: XQCoin^q = 0) = as the direct consequence of the Ward 
Identities discussed in section 2. 

The low energy description of the system in confined phase with a single quark is 
accomplished by introducing the r]' meson. As usual, the rj' would be conventional massless 
Goldstone boson if the chiral anomaly is ignored. In the dual sine-Gordon theory the rj' 
field appears exclusively in combination with the 9 parameter as 9 ^ 9 — r]'. As it is 
well known, this is the direct result of the transformation properties of the path integral 
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measure under the chiral transformations ijj — ?■ exp(i75y)'i/'. Therefore we have, 



N 
a=0 



(4.12) 



2L \2-kJ 



(Va) 



S^. = [ d'x.^ (Vr?') 



int 



N 



a=l 



COS I Oa ■ cr + 



e-7]' 

N 



where coefficient c determines the normahzation of the ry' field and has dimension one. This 
coefficient, in principle, can be computed in this model, but such a computation is beyond 
the scope of the present work. In four dimensional QCD the coefficient c is expressed in 
terms of standard notations as (c/L) — )• /^,. In terms of these parameters the rj' mass is 
given by 



cN' 



(4.13) 



Since ij' shows up in the Yang-Mills Lagrangian as r]'FF, we can compute our requisite 
expectation, {FF,FF), by functional differentiation with respect to rj' , 



(9(x)g(y)> 



1 i 5 i 6 
Z 5rj'{x) 5r]'{y) 



(4.14) 



Thus we have, 



(9(x)g(y)) 



1 



Z 6rj'{y) 
1 



Z 



VaVrj' 
(5(x - y 



VaVrj' 

c 



< 



N 

/ d^n 5{ri - x) 



(fri 5{yi - x)J(ri - y) 



N 



f r ^ 

/ d-^n / d^r2 6{ri - x)5{r2 - y) ^ 



a,b=l 



7]'{ri) r]'{r2) 

N N 



|(V(x)^'(y)) 



(4.15) 



The first term in (4.15) is precisely non-dispersive contact term with the "wrong sign" that 
we computed previously in pure gauge theory using two different methods, see (4.10) and 
(4.11). The second term represents the conventional dispersive contribution of the physical 
r]' state ^. One can compute it by redefining t]' — )• vj j ^fc field to bring its kinetic term 5^/ 



^This additional interactions due to the 77' exchange may in fact be used as a probe to study the relevant 
topological charges present in the system. It was precisely the idea behind the proposal, see relatively recent 
papers [15, 16] and earlier references therein, that Q jN behaviour unambiguously implies that the relevant 
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to the canonical form. In the lowest order approximation it is reduced to the conventional 
Green's function of the free massive rl scalar field with mass determined by (4.13), such 
that 



XQCD = I d'^x{q{x)q{y)) = ^ I d^^^ 



(5(x) — jn. 



2 



0, (4.16) 



where we represented the canonical -q' propagator in terms of its free Green's function in 
three dimensions. 

The structure of this equation follows precisely the same pattern we observed in analy- 
sis of two dimensional QED, see (2.6). Indeed, it contains the non-dispersive term due to the 
degeneracy of the topological sectors of the theory. This contact term (which is not related 
to any physical propagating degrees of freedom) has been computed using monopoles de- 
scribing the transitions between these topological sectors (4.10). The second term emerges 
as a result of insertion of the massless quark into the system. It enters xqcd precisely in 
such a way that the Ward Identity XQCoi'niq = 0) = is automatically satisfied as a result 
of cancellation between the two terms in close analogy with two dimensional case (2.9). We 
should also mention that very similar structure emerges in real strongly coupled QCD in 
the framework where the contact term is saturated by the Veneziano ghost. This structure 
has been confirmed by the QCD lattice studies, see [17] for the details and references on 
original lattice results. 

4.3 Interpretation 

The results derived in previous sections were formulated in Euclidean space using conven- 
tional Euclidean path integral approach. Our goal here is to give a physical interpretation 
of these results in physical terms of formulated in Minkowski space time. First of all, the 
(5(x) function which appear in the expression for topological susceptibility (4.10) should, 
in fact, be understood as total divergence, 



X 



I Si.) A = / A 8, = £ dE„ (j^^) . (4.17) 



Indeed, the starting point to derive x was the topological density operator (4.7) which is 
expressed in terms of 5(x — Xj) functions, but in fact represents the topologically nontrivial 
boundary conditions determined by the behaviour at a distant surface 52 as (4.17) states. 
The representation (4.17) explicitly shows that we are not dealing with ultraviolet (UV) 
properties of the problem wherein our approximation breaks down. Our treatment of the 



vacuum fluctuations must have fractional topological charges 1/A''. In the present weakly coupled regime 
these ideas have a precise realization as the basic vacuum fluctuations are indeed the fractionally charged 
monopoles (4.7). The results of [15, 16] are in fact much more generic as they are not based on a weakness 
of the interaction or semiclassical expansion, but rather, on generic features of the 77' system which are 
unambiguously fixed by WI. In the approach advocated in [15, 16] one can not study the dynamics of 
fractionally charged constituents in contrast with the present paper where the dynamics is completely fixed 
and governed by (4.12). However, the fact that the constituents carry fractional topological charge 
can be recovered in the approach [15, 16] because the color-singlet 77' field enters the effective Lagrangian 
in a unique way and serves as a perfect probe of the relevant topological charges of the constituents in the 
system. 
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problem is perfectly justified as (5(x — Xj) functions actually represent the far infrared (IR) 
part of physics rather than UV physics. This explains why our description is valid for 
X ^ L in spite of presence of the apparently UV singular elements such as the (5(x — Xi) 
functions which appear in equations of sections 4.1, 4.2. 

Our next comment is about the interpretation of the classical monopole gas from sec- 
tions 3. The monopoles in our framework are not real particles, they are pseudo-particles 
which live in Euclidean space and describe the physical tunnelling processes between differ- 
ent winding states \n) and |n + 1). The grand canonical partition function written in terms 
of the classical Coulomb gas (3.33) is simply a convenient way to describe this physics 
of tunnelling. In particular, the monopole fugacity (" together with factor should be 
understood as number of tunnelling events per unit time per unit volume 



where extra factor N in (4.18) accounts for N different types of monopoles present in the 
system. The expression (4.18) is precisely the contact term, up to factor computed 
in (4.10). It is not a coincidence that number of tunnelling events per unit time per unit 
volume precisely concurs with the absolute value of the energy density of the system (4.11), 
since the energy density (4.11) in our model is saturated by the topological fluctuations 
which are not related to any physical propagating degrees of freedom. We emphasize that 
while this energy density is not related to any fluctuations of real physical particles, this 
energy, nevertheless, is still real physical observable parameter, though it can not be defined 
in terms of conventional Dayson T -product. Instead, it is defined in terms of the Wick's 
T product, see Appendix of ref.[8] on a number of subtleties with definition of the energy. 

Finally, the characteristic Debye screening length which appears in the Coulomb gas 
representation in section 3 



should be interpreted as a typical distance in physical 3d space in which the tunnelling event 
is felt by other fields present in the system. The tunnelling interpretation also explains the 
"wrong" sign in residues of the correlation function (4.10) as we describe the tunnelling 
in terms of the Euclidean objects interpolating between physically equivalent topological 
sectors |n) rather than the tunnelling of conventional physical degree of freedom between 
distinct vacuum states in condensed matter physics. 

5. Conclusion and future directions 

The main results of this work can be formulated as follows. We studied a number of dif- 
ferent ingredients related to 6 dependence, the non- dispersive contribution in topological 
susceptibility with the "wrong sign" , topological sectors in gauge theories, and related sub- 
jects using a simple "deformed QCD" . This model is a weakly coupled gauge theory, which 
however has all the relevant essential elements allowing us to study difficult and nontrivial 




(4.18) 




g 



> L 



(4.19) 
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questions which are known to be present in real strongly coupled QCD. Essentially we 
tested the ideas related to the U{1)a problem formulated long ago in [1, 2, 3, 4, 5, 6] in 
a theoretically controllable manner using the deformed QCD as a toy model. One can 
explicitly see microscopically how all the crucial elements work. 

As this model is a weakly coupled gauge theory, one can try to formulate (and an- 
swer) many other questions which are normally the prerogative of numerical Monte Carlo 
simulations. One such question is the study of scaling properties of the contact term. We 
can address what happens to the contact term when the Minkowski space-time R^^i gets 
slightly deformed. For example, what happens when infinite Minkowski space-time R^^i is 
replaced by a large, but finite size torus? Or, what happens when the Minkowski space- 
time i?3^i is replaced by FRW metric characterized by the dimensional parameter R 
describing the size of horizon (H here is the Hubble constant)? A naive expectation based 
on common sense suggests that any physical observable in QCD must not be sensitive to 
very large distances ~ exp(— Aqcd-R) as QCD has a mass gap ~ Aqcd- Such a naive 
expectation formally follows from the dispersion relations similar to (2.2) which dictate 
that a sensitivity to very large distances must be exponentially suppressed when the mass 
gap is present in the system, and there are no any physical massless states in the spectrum. 
However, as we discussed in this paper, along with conventional dispersive contribution 
(2.2) in the system, there is also the non-dispersive contribution (2.3) which emerges as a 
result of topologically nontrivial sectors in four dimensional QCD. This contact term may 
lead to a power like corrections R~^ + 0{R)~^ rather than exponential like exp(— Aqcd-R) 
because the dispersion relations do not dictate the scaling properties of this term. In fact, 
this term in "deformed QCD" in infinite Minkowski space has been explicitly computed 
in this paper and it is given by (4.10), (4.16). As our model is a weakly coupled gauge 
theory, one can, in principle, compute the correction to the formulae (4.10), (4.16) due to 
the finite size of the system [18]. In different words, one can try to compute the corrections 
to the monopole fugacity C when the model is formulated in a finite manifold determined 
by size R. We suspect that the correction may demonstrate a Casimir power like behaviour 
R~^ + 0{R)~'^ rather than exponential like behaviour exp(— Aqcd^) similar to the previ- 
ously studied case of two dimensional QED [19] wherein all relevant ingredients are present 
in the system as discussed in section 2. Such a Casimir type behaviour may have profound 
consequences for cosmology as discussed in [8]. 
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